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Abstract 

The interpretation of the linear character of the observable classic rectilinear diameter law and 
the linear character of the Zeno-line (unit compressibility line Z = 1) on the basis of global 
isomorphism between Ising model (Lattice Gas) and simple fluid is proposed. The correct definition 
of the limiting nontrivial Zeno state is given and its relation with the locus of the critical point 
is derived within this approach. We show that the liquid-vapor part of the phase diagram of the 
molecular fluids can be described as the isomorphic image of the phase diagram of the Lattice Gas. 
It is shown how the the position of the critical points of the fluids of the Lenard-Jones type can be 
determined basing on the scaling symmetry. As a sequence the explanation of the well known fact 
about "global" cubic character of the coexistence curve of the molecular fluids is proposed. 

PACS numbers: 05.70.Jk, 64.60.Fr, 64.70.F 
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INTRODUCTION 



Since the seminal doctoral thesis of van der Waals [|TJ where the Principle of the Cor- 
responding States was formulated the search for the unifying principles for the description 
of the variety of thermodynamical properties of complex matter is the key point of statis- 
tical physics. By the rigorous methods of the Statistical Mechanics the universal character 
of the van der Waals (vdW) equation was demonstrated and the numerous extension were 
proposed 0. The consequences derived from the vdW equation and the observed devia- 
tions became the driving force for the further theoretical developments. One of the main 
achievement of van der Waals was the description of the phase liquid- gas equilibrium which 
terminated at the corresponding critical point (CP). It opened the possibility to connect the 
observable thermodynamic quantities with the characteristics of intermolecular interactions. 
In the modern theory of the Critical Phenomena the ideology of the isomorphism classes of 
the critical behavior provides the description of the real systems using the results obtained 
for the model systems like lattice models ||. In particular the molecular liquids with short 
range interactions of the Lenard- Jones (LJ) type belong to the isomorphism class of the 
Ising model. The latter can be easily mapped to the Lattice Gas (LG) model, which is 
determined by the Hamiltonian: 

H = - J^mnj - ft m , (1) 

< ij ) i 

where n, = 0, 1 whether the site is empty or occupied correspondingly. The quantity J is 
the energy of the site-site interaction of the nearest sites i and j, jj, is the chemical potential. 
The temperature variable corresponding to the Hamiltonian (|XJ) is denoted as t. The order 
parameter is the probability of occupation of the lattice site x — (rii) and serves as the 
analog of the density. The phase diagram of the LG is symmetrical with respect to the line 
Xq = 1/2 and formally extends up to the low temperature region t — y 0, where the limiting 
states x = and x = 1 exist only. 

Two of the consequences of the vdW equation are of primary importance for our consid- 
eration. The first one is the well-known rectilinear law for the diameter of the coexistence 
curve in terms of the particle density n and the thermodynamic temperature T: 



where rii, i = liq, gas are the values in the corresponding phases, n c ,T c are the critical 
density and the temperature. The second fact is another linear law which was derived by 
Batchinsky j3J long ago from the classical vdW equation: 

P = --an 2 . 3 

l-nb v ; 

The curve defined by the equation Z = P/(nT) = 1 for the compressibility factor Z as it 

directly follows from Eq. (Q) is a straight line: 

n T 

YJb + T £dw) = 1 ' ( > 

where j^ vdW ) = a/6 is the Boyle temperature determined by the vdW constants in accor- 
dance with the virial expansion for Eq. (JJ). In work of Ben-Amotz and Herschbach || the 
line Z — 1 is called by the Zeno-line. 

Both relations (Q) and (|5|) are approximate and have phenomenological character. But 
surprisingly they are observed for a wide variety of fluids. As was noted in || the law 
(0) is fulfilled in rather broad temperature interval of the vapor-liquid phase coexistence 
region 0.65 < T/T c < 1. In the vicinity of the CP the deviations from the linearity become 
noticeable [[7|, [|. 

In series of works @-13| the performed analysis of the data for both the real fluids and 



the systems with model potentials discovered a number of interesting relations between 
coordinates of the CP, the law of the rectilinear diameter @ and the Zeno-line: 

n T 

- + — = 5 

n b T B 

Here Tg is the Boyle temperature and for 7i& the following equation followed from the virial 
expansion was proposed in [[L4| : 

The phenomenological concept of the "Triangle of Liquid-Gas States" has been formulated 
in [11 1 . It is heavily based on the cute observation that the Zeno-line is the tangent to the 
binodal extrapolated to low temperature region. Then the liquid-vapor binodal is inscribed 
into triangle formed by the straight lines: the coordinate axes n and T and the Zeno-line. 

Obviously, the thermodynamic states which satisfy the equation Z = 1 include the trivial 
states with negligible density n — > 0. This is nothing but the coordinate axis on the ther- 
modynamic n — T plane. Thus the set of the Zeno states is formed by the union of two lines 



n = and (|5]). The statement ([5]) about the linear nature of Zeno-states is far from trivial. 
Indeed the expression for Z, which for the system with binary potential interaction $(ri2) 
reads as (see e.g. 0): 

Z = — = 1-— F r— g 2 {r;n,T)dr, 7 

nl 61 J or 

where g 2 is the pair correlation function. Therefore the linearity of the Zeno states is 

due to quite specific structure of the dependence of correlation function g 2 (r;T,n) on the 

thermodynamic state. 

Below, following the arguments of previous work Hl5l , we show that these facts can be 
casted into elegant geometrical formulation and expand them into general case of the short 
range power-like attractive potentials in d dimensions. We show that the locus of the CP can 
be estimated using the correspondence between the scaling properties of the Hamiltonians 
for the lattice gas and the fluids with the power-like interactive potentials. 



THE ZENO-LINE AND THE GLOBAL ISOMORPHISM 



As was shown in |15|] the linearities of (0) and (|^) can be derived on the basis of the 



assumption of the existence of the global isomorphism of the real liquid-gas part of the 
phase diagram of the lattice gas model. Note that the line x — 1 of the LG can be thought 
of as the analog of the Zeno-line. Indeed the pair correlation function g 2 (r) of the LG 
vanishes identically for such "holeless" states according to the definition. The same is true 
for the state with empty sites x — 0. Thus the line x = can be identified with the zeroth 
density axes n = of the real fluid, while the line x = 1, obviously, can be identified with the 
Zeno-line Z — 1. To conserve the linearity and the adjacency properties of the characteristic 
elements the isomorphism between LG and the fluid should be chosen in the class of the 
projective mappings [I6j and has the form: 

x m m zt 



where 



T 

(9) 



T z -T c 

The coordinates of the CP for the liquid are: 



T c = T z —— , (10) 



2(l + z)' l + z 



provided that the LG temperature variable t is scaled so that t c = 1. In view of Eq. (|9D it 
should be noted that if Tz and rib are fixed then the parameter z parameterizes the locus of 
the CP: 

n c T c _ z+1/2 

n b T z l + z 1 ' 

This means that for the substances belonging to the same class of the corresponding states 
the loci of the critical points scaled to the Tz and nj, lie along the straight line. This 



correlates with the empirical analysis in ||14|| . Because both Tz and T c are determined by 



the interparticle potential the parameter z can be connected with its symmetry properties. 
Note that all results obtained rely solely on the geometrical facts about phase diagrams of 
the Lattice Gas and the fluid and do not depend on the specific details of the interactions. 
Two limiting Zeno states T — > 0, n — > n& and T — > Tz, n — > form the triangle of the 



Zeno states |lT|. It should be noted that the temperature Tz of the limiting Zeno state 
with nz — > does not necessary coincide with the Boyle temperature T B as was originally 
proposed in ||. In || the Boyle state n — > , T = T B as an obvious limiting Zeno state is 
used and then the parameter rib is defined by Eq. (U). But any state on n-axes can be used 
in such way. Two limiting Zeno states n — > , T = Tz and n — n& , T — y belong to two 
physically different thermodynamic regions. Therefore there is no apriori physical reasons 
to define the value rib basing on Tg. From the point of view of the virial expansion: 

oo 

z = ^f = l + Y. B ^ T ) nk ( 12 ) 
k=i 



where Bk(T) is the virial coefficients [17[, the limiting state Tz with nz — > is the point 
where the nontrivial branch of the solution of Z = 1 emerges from the axis n = with the 
specific tangent determined by T z itself. Clearly the Boyle temperature T B determined by 
the condition: 

B 2 (T B ) = 0, (13) 
is the simplest choice but in general T z ^ T B because flT3"|) does not guarantee the correct 



slope of the Zeno line. Indeed if one search the solution of the equation: 

Z-l=nB 2 (T) + B 3 (T)n 2 + ... = (14) 

in a form of 



in the vicinity of T = T z then obviously the infinitesimal condition: 



dZ{n z {T),T) dZ 



dT dn 



dn z d Z 
+ 



T->T z ,n^0 dT dT 



= B' 2 (T) + B 3 ^ = (16) 



should be fulfilled. Using Eq. fll2f) we get that T z is the root of the equation: 

Tz= " b lm)' (17) 

This is the rigorous definition of the limiting Zeno state at n — > from which the Zeno line 
emerges. Note that for vdW equation with B 2 = b — f;, B 3 = b 2 and n b = 1/6, Eq. flT?D 
gives the standard result T z = a/b. Therefore this is the specific property of the vdW 
equation that the temperatures T z and Tb coincide. Actually the value T z depends on the 
dimensionless packing parameter rjb = rib & d i n d dimensions, where a is the characteristic 
scale of the potential, e.g. the hard core diameter. This parameter is determined for the 
opposite limiting Zeno state with T — > 0. Though Eq. (|T7p formally coincide with Eq. (^), 
the approach proposed for the determination of the limiting Zeno-state at n — > is different 
from that used in f9j. Thus the situation is inverse: rather the value n& determines the value 
T z and not vice versa. In order to make this consistently with the definition of the Zeno-line 
we choose so that T z = T B vdW \ where T^ dW ^ = a j]j j s the Boyle temperature in the van 
der Waals approximation. We use this to calculate the locus of the critical point in Section . 
Using the standard expressions for the virial coefficients B^ for Lenard- Jones potential 

m= - Wo ((°-y-(°-yy (18) 

the solution of (|TT| ) for 3-dimensional system gives a 3 0.965 and Tz = 4 while Tb ~ 3.42. 



SCALE INVARIANT MEAN-FIELD THEORY 



In this section using the proposed isomorphism we propose the variant of the mean-field 
approach to calculate the locus of the critical point. It exploits explicitly the scale symmetry 
which inherent to the power-like interaction potentials. We call it scale-invariant mean-field 
theory. 

Let us show that the parameter z can be determined if the attractive part of the inter- 
action potential possesses the scaling symmetry. The starting point is that using ( |T0"1) it is 
easy to derive the following relation: 

d In ( TJT Z ) = 1 
d In ( n c /rib ) z 
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The locus of the CP is determined mostly by the attractive part of the interaction potential 
and the size of the particles (see e.g. flTf ). These parameters can be connected with the 
corresponding ones for the LG model. Suppose that the potential of the <i-dimensional 
system is $(r) = $ rcp (^) + $ a ttr( r ), where $ rep , $ a ttr are the repulsive and the attractive 
part correspondingly. Additionally we assume that the attractive part has the power-like 
behavior $ rep (r) ~ r -(d+e) £ > q Then the corresponding energy of the interaction for the 
configuration of the number density n(r): 

E int = y / $ ( ri2 ) n ^ n{r ^ dr ™ ( 2 °) 

can be compared with the LG Hamiltonian (P. Now we put the constraint of the scale 
invariance of the partition function: 

N 

Z= e~ E ^' T dT, dT = Y[ dn , (21) 

"* i=l 

corresponding to ( f20|) with respect to scaling transformations of n — > X An n and T — > X At T, 
where A is the scale parameter. The exponents A n and can be found using standard 
similarity considerations widely used in condensed matter theory (see e.g. [0). E.g. for 
the coordinates of the CP we have: 

1 



//,. ~ - T c ~ $ attr (r c ) = $ attr {n- 1 ' d ) , (22) 



-l/d 

where r c is the mean interparticle distance in the critical state. Therefore from (|T9|) and 



c 



([22|) we can write the relation: 

d In ( TJT Z ) d In $ at tr(n c ~ 1/rf ) (23) 



d\nn c /rib d\nn c /nb 
which gives the following consistency relation: 

l/z = l + e/d. (24) 

Thus we obtain the one-parameter group of scaling transformations for the locus of the 
CP connected with the change of the characteristic scale of the interaction: 

n c {z) -> n c (0) A^ 1 , T c (z) T c (0) A( 1+£ / d ) . (25) 

Eq. ( p5j) is equivalent to the constraint: 

n c {\) 1+£ / d T c {\) = const. (26) 



Let us show that Eq. ( p4[) can be considered as some consistency condition between the 
scaling of for the LG Hamiltonian (|1|) and that for (^). Note that Eq. fl25|) represents in 
particular the scaling of the molecular size n b and corresponding scaling of energy scale given 
by Tg or equivalently Tz- To connect the LG model with the LJ fluid where the interaction 
has the form ( |2"0"| ) we take into account that under the scale transformation the loci of 
the corresponding critical points of these systems should be scaled consistently. Indeed, 
according to the structure of the LG Hamiltonian (0) the analog of Eq. (|2~6|) for the critical 
point of the LG is as following: 

x\ t c = const — ^ = -2 — . (27) 

Cii X Q X Q 

Then the consistency condition between the scaling of the critical point coordinates of the 



LJ fluid flTJ) and the Lattice Gas (|27|) is: 



TJT *'l + t)=2 * !=! + £ (28) 



(29) 



n c /rib V d ) z d 

In accordance with (§) we get: 

T c _ 1 n c 1 + | 



T Z 2 + §' n b 2(2 + |)- 

The case of the attractive potential of the van der Waals forces in <i-dimensions corresponds 
toe = d,z = l/2 and we get the estimates: 

T c = T z /3 , n c = n b /3 . (30) 

In the following Section we use ( p9|) to calculate the locus of the critical point of the fluid 
where the interaction potential is of the LJ type (|18|). 



THE LOCUS OF THE CRITICAL POINT FOR THE LENARD-JONES FLUID 



To test the predictions (|29|) we compare them with the results of the computer simulations 
for the potential (|18D for 2D, AD and 5D dimension geometries available in [|19[ [20|1 . 



In accordance with the results of Section we determine the packing parameter n b = n b o~ d 
so that: 

Tz = T^ = a/b=-^-, « = 2^^, b = 2 d -\ (31) 

o — a o — a 

where the particle volume Vd(a) and U are put to units. The result are placed in Table |. 
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LJ "6-12" fluid 


2D 


3D 


4D 


5D 


T z 


2 


4 


8 


20 


rib & d 


0.91 


0.965 


1.342 


2.54 


T B 


1.56 


3.418 


9.01 


40.4 



TABLE I: The results of computation of Tz and Tb for LJ potential ( |l8| ) for different dimensions 
and the packing factor re& a d . 

The locus of the CP is obtained easily according to (p9|) and (|3lD (see Fig. [DID- The 
comparison with the results of the simulations is in Table 0. Note that expressions (^) and 



( j31~l) allow to clarify the fact noted in |20| about different characters of the dependencies of 
T c and n c on the dimensionality d. The value T c strongly depends on the dimension because 
of the dependence of T B while the ^-dependence of n c is rather weak. 

Note that nonmonotonic dependence of the critical density on the number of dimensions 
d can be interpreted in terms of the proposed isomorphism as follows. The equilibrium 
interparticle distance for the LJ potential (|T8|) is r = <r2 1//6 and can be considered as the 
spacing of the cubic lattice for the LG. The effective radius r e g of excluded volume occupied 
by the particle from the physical point of view can be defined by the obvious energetic 
condition: 

$(2r eff )=T r cff (T). (32) 

So the packing of spheres of radii r eff on the corresponding cubic lattice can be found as 



n(T,d) 



2 d/6 a d 



V{d)r^{TY' 



(33) 



where V(d) = 2^^ is the volume of the unit ball in d dimensions. Obviously n c depends on 

1 ("2 ) 

rj(T, d) monotonically. The dependence of r](T, d) on the dimension d is shown on Fig. |2] by 
grey lines and shows the minimum in the interval d < 4 in dependence on the temperature. 
Note that at T = 0, when r e fj = a this minimum reaches exactly at d = 4. 

The existence of global isomorphism also provides the explanation of the fact about the 
cubic form of the binodal for real molecular fluids. The fact that the shape of the binodal 
is almost cubic and characterized by the effective exponent /3 e // = 1 ~ 1/3, ip = n '~ ng 
in a wide temperature interval for the system with short range interactions was discovered 

(see also 



long ago by J. E. Verschaffelt [21H 



and since then has been confirmed by 
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LJ "6-12" fluid 


2D 


3D 


4D 


5D 




0.5 


1.333 


3.2 


9.08 


r c (nnm) , g 


0.515 


1.312 


3.404 


8.8 (?) 


n c 


0.353 


0.322 


0.404 


0.693 


n { r m \ |2C] 


0.355 


0.316 


0.34 





TABLE II: The comparison of the estimates (^) with the results of numerical simulations [20]. The 
value T C 5D ps 8.8 of the critical temperature for 5D case was obtained in |^] by simple extrapolation 
of the results for lower dimensions. 
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FIG. 1: Critical temperature T c for the potential (jig) of the d-dimensional LJ fluid according to 
Eq. (pi). 



the direct numerical simulations |24] , Po| . In particular in the principle possibility to 
derive the nonclassical exponents for molecular liquids basing on the specific closure of the 
Ornstein-Zernike equation. 

The proposed global isomorphism between LJ fluids and LG leads to the conclusion that 
the phase coexistence region of the latter is nothing but isomorphic image of the corre- 
sponding region of the phase diagram of the LG. For example applying the isomorphism 
transformation (|8]) to the known result for 2D Ising model [p7[| : 

1 



x 



1/2 ± fit) 1 '*, /(f) = 1- 



sinh 4 ( 2J/t 



(34) 



10 



0.8 



0.6 



0.4 



0.2 - 



T=10 




FIG. 2: Critical density n c for the potential (18) of the d-dimensional according to Eq. (|29| ) (dashed 
line is the guide for an eye) and the packing fraction rj(T, d) given by Eq. (|33|) (grey lines) . 



T/T n 




FIG. 3: Binodal for 2D Lenard- Jones fluid with $ a ttr( r ) ~ r calculated as the isomorphic image 
of the binodal for 2D Lattice Gas (Ising model) in accordance with (|8]). 

the binodal for 2D LJ fluid can be easily obtained (see Fig. |3|) and compared with the results 
of the computer simulations |[28|| . Also the inference that the cubic character of the binodal 
for real systems is the consequence of the cubic character of the binodal for the Ising model. 
So one can expect that the coexistence curve of the 3D Ising model is described in the 
"zeroth approximation" by the algebraic cubic curve in analogy with (134]): 

(35) 



x=l/2±/(t) 1 / 3 ^t(x) = r 1 (|2x-l| 3 ) , /(1) = 



where f(t) is analytic function of the temperature. 

In particular the isomorphism transformation (0) applied to (P3|) allows to connect the 
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critical amplitude B for the order parameter of the LJ liquid: 

n ~ n " ±B^\rf + ..., (36) 



with that of the Ising model: 

2x-l = ±B$ s) \rf + ... . (37) 

Taking into account the relation (||) between temperature variables T and t via compar- 
ison (jSHj) and (|37D we obtain the following relation between amplitudes 

As the amplitude Bq Is ^ is known form exact result at D = 2 or from computer simulation 
for D = 3 (see e.g. [[^|) one get the critical amplitude for LJ fluid. For the real 3D fluids 
this gives Bq LJ ^ 1.69 for 2D case and Bq LJ ^ m 1.85 for 3D case. Then the proximity of 
the numerical values of the critical exponents for the systems of the Ising model universality 
class to the rational numbers is quite natural. They do not necessary close to the mean- 
field values of the Landau theory and for the systems with the short-range interactions the 
mean- field behavior is not observed far away from the CP |3"0||. This gives the grounds 



for the universal global cubic law for the binodal proposed in fbl. In essential, the same 



global cubic law is widely used in computer simulations of the system with short ranged 



interactions |29], pl| . The fluctuations renormalize these rational exponents near the CP 
through small exponents r\ and a, which determine the behavior of the correlation functions 
of the density and the entropy in the immediate vicinity of the CP [|3"2^ . The fluctuations 
are responsible for the deviation from the linear diameter law (fj). If we assume that the 
transformations between thermodynamic averages like x and n are generated by the trans- 
formations between the corresponding microscopic fields then such deviation is interpreted 
as the consequence of the nonlinearity of the transformation for the field variables. Indeed 
for the nonlinear functional F of the thermodynamic field n(r) in the presence of strong 
fluctuations (F[n(r)]) ^ F((n(r))), where (...) denotes the statistical average. From 
this point of view the transformation (|8|) can be considered as the "mean- field" analog of 
commonly used linear mixing of the fields in the vicinity of the CP |3^ . 
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CONCLUSIONS 



In concluding Section we summarize the main results and discuss possible restrictions of 
the proposed approach. 

The whole set of phenomenological facts about the interrelations between Zeno line, 
the rectilinear diameter, the binodal and the locus of the critical point revealed in [|10|- 
[L3| needs in unifying view from the microscopic approach of Statistical Mechanics. The 
proposed isomorphism allows to search such view in terms of the correspondence between 
the Hamiltonians of the LG and the fluids with short-ranged pair potentials. In particular 
it seems natural to connect the thermodynamic potentials for these systems. This allows to 
determine the parameters of mixing for the fields of the order parameter and the temperature 



33| . The results obtained for the locus of the critical point of LJ fluid and critical amplitudes 



of the equation of state show good correspondence with the results of computer simulations. 
Though the attractive part of the potential was taken into account in proposed version of the 
scale invariant mean-field approach the relations Eq. fl22|) , (|23|) can be generalized to include 
the repulsive part of the potential. E.g. instead of ([23|) the following relation arise: 

d In (T c /T z ) _ dln<l> attt (n- 1/d ) 
d In n c /rib d\nn c jrib 

This relation along with the first equation in ( |10D would give the closed system of equations 
to determine both n c and z consistently. For the considered case of the LJ potential the 
repulsive part gives negligible correction so we neglect it here. But in general case the 
influence of the repulsive part on the locus of the critical point should be taken into account. 

Note that the isomorphism relies on simple geometrical facts among which the tangency 
of the Zeno-line to the extrapolated liquid- vapor binodal into hypothetic region T — > is the 
most crucial |T3J . For real systems the liquid branch is restricted by the triple point and this 
causes the problem of the correct determination of the parameter rib- Obviously the slope of 
the Zeno line at rib influence essentially the position of the opposite Zeno state with n — > 0. 
To overcome this difficulty the correct determination of the limiting Zeno state (Tz, 0) was 
proposed along with the determination of rib form the constraint Tz = a/b where a and b 
are the corresponding van der Waals constants for the equation of states. In principle it is 
possible to determine the values of rib using the data for the locus of the triple point and 
the slope of the binodal using the isomorphism transformation which allows to relate this 
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data with the ones for the LG binodal. This will be the subject for further study. 

It is interesting to note that in the limit e — » — d, i.e. when the range of the interaction 
becomes infinite, T c — > Tb- Thus the coexistence curve degenerates into triangle. This 
situation can be compared with the phase diagram for the infinite diluted polymer with the 
Flory theta-point (see also in 
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